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Jorge Moreno 1 



O 

(N 

W) 
< 

o 

u 

43 
6 

(N 
> 

m 
m 
m 

O 



1 SISSA, Astrophysics Sector, via Bonomea 265, 34136, Trieste, Italy 



ABSTRACT 



Accumulating evidence suggests that galaxy interactions play an important role in 
shaping the properties of galaxies. For this reason, cosmological studies focused on the 
evolution of halo/subhalo pairs are vital. In this paper I describe a lar ge catalogue of 
halo pairs extracted from the publicly available Millennium Simulation (ISpringel et al.l 
l2005f ) . the largest of its kind to date. (Throughout this work I use the term 'halo' to 
refer both to individual haloes in the field and subhaloes embedded in larger struc- 
tures.) Pairs are selected according to whether or not they come within a given critical 
(comoving) distance d C rit, without the prerequisite that they must merge. Also, a con- 
dition requiring haloes to surpass a critical mass M cr - lt during their history is imposed. 
The primary catalogue, consisting of 502,705 pairs, is selected by setting d cr it = 1 Mpc 
hr 1 and M cr it = 8.6 x 10 10 Af Q /j -1 (equivalent to 100 simulation particles). One of the 
central goals of this paper is to evaluate the effects of modifying these criteria. For this 
purpose, additional sub-catalogues with more stringent proximity and mass conditions 
are constructed (i.e., d cri t = 200 kpc ft," 1 or/and M crlt = 8.6 x 10 11 M s ft,- 1 = 1,000 
simulation particles - see Table [T] for a summary). I use a simple five-stage picture 
to perform statistical analyses of their separations, redshifts, masses, mass ratios, and 
relevant lifetimes. The fraction of pairs that never merge (because one of the members 
in the pair is absorbed by an third halo or both members survive until the present 
time) is accounted for. These results provide a broad picture that captures the es- 
sential characteristics behind the evolution of these halo pairs. This is the first of a 
series of papers aimed to explore the huge wealth of information encoded in this cata- 
logue. Such investigations will play a fundamental role in future cosmological studies 
of interacting galaxies and binary (and multiple) quasars. 

Key words: cosmology: theory - dark matter - large-scale structure of the Universe 
- galaxies: formation - evolution - interactions - statistics 



1 INTRODUCTION 



Galaxy mergers and interactions are among the most spec- 
tacular events in the sky. But more importantly, they pro- 
vide a bridge between worlds: linking the large-scale hierar- 
chical universe to the internal properties of individual galax- 
ies and their evolution. Such encounte rs have the ability 
to ch ange the morphology of a galaxy (|Toomre fc Toomrd 
1 19721 ), and can cause a substant ial migration of its gas to- 
ward s the central regions (e.g., iBarnes fc Hernquistl Il99ll . 
1996). Ultimately, sudden bursts of star formation and the 
activation of the supermassive black hole at the nucleus 
can be attributed to this mechanism. This single idea has 
been, without a doubt, instrumental to all modern theories 
of galaxy formation (see e.g, the textbook bv lMo et al.ll201ol . 
and references therein). 

Admittedly, mergers have received more attention than 
the early stages of interaction preceding them. The last 



few years have witnessed the development of a whole in- 
dustr y of pre-prepared galaxy-galaxy merger si mulations 
(e.g., iHopkins et all 120051 : iRobertson et al.l |2006| U Unfor- 
tunately, the pre-merger stages are sel dom the prior- 
ity o f these numerical studies (b ut see |Pi Matteo et al.l 
120071 ; ICallegari et al. I 120091 . |2011| . whose simulations are 
tailored to emphasize the role of early interactions on 
star formation and black hole growth respectively). Like- 
wise, huge efforts have also been reserved towards cos- 
mological theoretical merger ra tes of galaxies , dark mat 



ter haloes and subhal o es (e.g., lLacev fc CoL 




fc 1993, 



1994; 



Governato et all Il999l: iGottlober et all l200lT ICohn et ail 



Benson! 120081: iGuo fc Whitd l200Sl: iFakhouri fc M, 



Wetzel et aJJ 



Neistein fc Dekell |20Qg|; iMateuj I200I" . _ 
Genel et alj|2009l: lMorendl2010l : Hester fc Tasitsiomil 
Hopkins et al]|2010f ). 



Recently, however, a wide array of observations have 
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provided hints that persuade us to redirect some of our 
attention to galaxy interactions. For example, a number 
of studies demonstrate that star formation is enhanced in 
pairs (typically within < 50 kpc), an d this effect is anti- 
correlated with their separation (e.g., iBarton et al.l 2000: 
Lambas et aDl2003l: iNikolic et al] 120041 : lEUison et aljbooa 



Pare et all l20ld : I Wong et all 1201 ll . and subsequent works 



by these groups). Even at high redshifts, Lya emission is 
heightened in pair s of Lyman-brea k galaxies (LBGs) with 
small separations (ICooke et a.l.ll201oh ■ 

Galaxy interactions als o play a role i n act ive galactic 
nuclei (AGN). For example. lAlonso et alJ l|200l ) report en- 
hanced nuclear activity in ~10% of their galaxy pair cata- 
logue (with projected separations < 25 kpc /i -1 ). The in- 
tensity of nuclear activity may depend on various factors, 
including the budget of gas available and the strength of 
the torques that can drive this gas into the central regions 
(the latter effect is primarily determined by the orbit of the 
encounter and the mass ratio). Interestingly enough, they 
find that galaxies with bright /massive companions appear 
to be more active, suggesting that galaxies with shallower 
potentials are more susceptible to the gravitational torques 
produced during the interaction. 

Lately, the search of binary quasars has gained impe- 
tus greatly in part because these systems are thought to 
be the precursors of coalescing black holes. While it is true 
that supermassive black hole growth is expected to be domi- 
nated by gas accret ion during quasar phases, not by mergers 
of black hole pairs JSalucci et al J 1 19991 : iMarconi etalll2004l ; 
IShankar et alJ l 20041200911 : these unique events remain inter- 
esting, especially because their existence could be detected 
by future gravity- wave missions like the Laser Inte rferom- 
eter Space Antenna (LISA) (e.g.. ISesana et al"1l2005l ). From 
a theoretical perspective, pairs of quasars at galactic sep- 
arations establish the initial conditions necessary to com- 
prehend these peculiar duos as they find their way into the 
central regions and perhaps into each other. 

In particular, b uilding up on p reviou s samples 
(|Hennawi et all 120061 : iMvers et all 12001 120081 ') . Foreman, 
Volonteri & Dotti (2009) compiled a list of 85 quasar bi- 
nary pairs. Their physical-separation distribution peaks at 
~30 kpc, with a long tail at larger separations (~ 100-200 
kpc). Indeed, these authors find 14 pairs separated by dis- 
tances > 300 kpc! The size of this cens us was recently in - 
creased by over an order of magnitude by I Liu et all (|2011bT ). 
who report 1,288 binary AGNs (256 with tidal features), 39 
triplets, 2 quadruplets and 1 quintuplet. For the binaries, 
their projected-separation distribution exhibits peaks at ~5, 
20 and 60 kpc h~ l . However, if only pairs displaying tidal 
features are considered, a sharper peak at ~20 kpc ft -1 re- 
mains (at this stage, these are the only ones guaranteed to be 
truly interacting, and further investigations might be desir- 
able for the rest). Additionally, these authors find evidence 
that the galaxies in the pairs with the smaller stellar mass 
manifest more nuclear activity than their companions, sug- 
gesting that either secondary galaxies have more available 
fuel or that they are more susceptible to tidal perturbations 
|Liu et al.ll201ld ) . Lastly, supermassive black hole pairs have 
also been observed in dust-en shrouded merging galaxies. For 
instance. iDasvra et all (|2006h find 21 ultraluminous infrared 
galaxies (ULIRGs) with nuclear separations between 1.6 and 



23.3 kpc, thought to be between the first and final passage 
prior to a merger. 

All this emerging observational evidence highlights the 
immediate need to focus on the long-lived evolution of 
galaxy pairs (and their supermassive black holes) in the 
stages prior to mergers. With this ultimate goal in mind, 
I analyse a very large catalogue of halo pairs and their evo- 
lution, suitable for future statistical studies of galaxy and 
quasar binaries. This catalogue, consisting of over half a 
million pairs, is extracted fr om the publicly avai lable Millen- 
nium Simulation Database |Springel et alj|2005f ). the largest 
of its kind published so far. In this work, pairs of haloes 
are selected by whether or not they are identified within 
a critical separation at some stage of their history. With 
this choice, pairs of haloes that never merge are not au- 
tomatically rejected by construction. The effects of chang- 
ing this critical condition are considered. I also explore the 
consequences of selecting pairs involving haloes that reach 
masses above a given threshold at some point of their evo- 
lution. The motivation here is to test whether close pairs of 
massive haloes are more or less susceptible to merging than 
pairs involving haloes with small masses. 

This is the first of a series of papers that exploit this cat- 
alogue with the aim to understand the evolution of halo pairs 
in the universe. As in any exploratory paper, I only touch 
on the most basic general trends. These results will serve as 
a benchmark for future studies - with more specialised sub- 
catalogues, refined recipes, and astrophysical applications 
- to be presented in forthcoming papers. Broadly speak- 
ing, I describe the evolution pairs in terms of five simple 
stages: initial, entry, closest, final and fate (see Section [2.31 
below for definitions). Quantities like the relative distance, 
redshift, total mass, mass ratio and lifetime are analysed 
statistically. Kinematic vector quantities like velocities and 
spins will be analysed in a companion paper (Moreno 2011b, 
in preparation). Relevant time intervals are also studied. Of- 
ten throughout the paper I use the term 'interacting' to refer 
to pairs within a chosen critical separation. It is important 
to emphasize that I do not presume to know a priori the 
exact proximity criteria at which tidal disturbances are ac- 
tually effective. Nevertheless, using this term facilitates the 
discussion of the trends in the halo-pair catalogue. 

The remainder of this paper is organised as follows. Sec- 
tion [2] describes the methods: the Millennium Simulation, 
halo-pair identification, the catalogue and the stages of evo- 
lution. Statistical results are found in Section [3] which fo- 
cus the stages and relevant lifetimes. In Section 2] I discuss 
the results, their interpretation, and their relation to other 
works in the literature. A summary is provided in Section [5] 
Throughout this work, all distances are comoving, unless 
stated otherwise. 



2 METHODS 

2.1 The Millennium Simulation 

In a nutshell, the Millennium Simulation is a dark- 
matter only simulation that adopts a fiat ACDM cosmol- 
ogy with the following parameters: (O m , ^b, o"8, n B , h) = 
(0.25,0.045,0.9,1,0.73). The simulation covers a periodic 
comoving box with Lbox = 500 Mpc h^ 1 on the side. It 
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follows the evolution of 2160 3 ~ 10 10 dark matter par- 
ticles with mass m p = 8.6 x W 8 M t r ) h~ 1 . Dark matter 
haloes are identified usi ng a friends-of-friends (fof) algo- 
rithm (|Davis et alj|l985h with a linking length set equal to 
0.2 times the mean inter-particle separation. These struc- 
tures are identified at 64 snapshots from z = 127 to the 
present, spaced nearly equally in log(l + z). The Plummer- 
equivalent force softening scale is set to 5 kpc h~ x . 

Self-bound structures with the FOF groups ar e iden - 
tified with the subfind algorithm (|Springel et al.l l200l[ ). 
This method allows to identify the main background dark 
matter halo a nd it s substructures (i.e., the subhaloes). 
iGiocoli et all (|2010h compares subfind against various 
other substructure-finding algorithms in their Appendix A. 
Merger history trees are built from these self-bound ob- 
jects by tracking their 'descendants' at the next snap- 
shot wi th lower redsh i ft (se e the supplementary informa- 
tion in ISpringel et al. 1 120051 . for more details). Once the 
trees are built, one can track the identity of each object 
in the hierarchy back in time by following its principal 
branch (see Figure [T] and its associated discussion below). 
This is roughly done by listing the most-massive progeni- 
tor, its_sjibsec[uent_jnost-ma^ive progenitor, and so forth. 
See [51 Lucia fe Blaizotl (|2007l ) for a discussion of the sub- 
tleties in assigning the most-massive progenitor to an object 
and their improved 'first progenitor' designation (hereafter 
I drop the word 'principal' unless it is absolutely necessary). 

A few comments on terminology are in order. In this 
paper, the words 'subhalo' or 'satellite' always refer to a self- 
bound object immersed in another larger self-bound struc- 
ture (the latter is often called the 'background halo' or the 
'host halo'). Objects that are not substructures of other ob- 
jects are called 'isolated' haloes. For instances, subhaloes 
were isolated haloes before becoming members of a larger 
halo. By definition, host haloes are isolated haloes with sub- 
structure. In this work the term 'halo' refers to any self- 
bound structures identified by the the subfind algorithm, 
without regard to whether they live in a host background 
halo or not. An upcoming paper will exploit these environ- 
mental differences. 

Lastly, defining mergers is not entirely trivial in N-body 
simulations because a subhalo may go below resolution be- 
fore having a chance to merge with another structure. Many 
authors adopt analyti c dynamical-friction recipe s to track 
down mergers (see e.g- lKitzbichler fc Whitel2008l . and refer- 
ences therein) . For the sake of simplicity, if two haloes share 
a descendant at the next snapshot, they are considered as 
merging. More sophisticated dynamical-friction techniques 
can always be incorporated in future applications of this 
catalogue. See Section [3] for more comments on this issue. 

2.2 Identifying the Halo Pairs 

The goal of this Section is to establish the language neces- 
sary to describe the evolution of close pairs of haloes. To do 
this, I discuss a single merger history tree extracted from 
the Millennium Simulation in detail. In reality, the spatial 
evolution of a merger tree should be represented in (3 + 1)- 
dimensions. However, for illustrative purposes, the (1 + 1)- 
dimensional representation in Figure Q] suffices. 

In this Figure, the black open circles depict the red- 
shifts and positions (schematically represented by the X- 



6 



4 



2 







1 1 1 1 1 1 1 1 1 1 
1 1 

d crit ["" 

i 

i 


1 1 
A B 

..<M>... 

1 I 

w 


1 I I I 1 I 

D 

C / 


- 
- 
- 




k 






• Merger 
Interaction 


V 






— Branch A 








. Branch B 




%T E 




Branch C 

Branch D 

Branch E 

, , , 1 , , , 1 , , 


1 I 


, I ?Sk« I , , 





0.2 0.4 0.6 0.8 1 1.2 
X [Mpc h" 1 ] 



Figure 1. Interactions in a merger history tree: redshift versus lo- 
cation (depicted schematically by the X-coordinate). Haloes iden- 
tified in the simulation are indicated by the open black circles, 
with radii proportional to the logarithm of the mass. Mergers are 
represented by solid red circles. Interactions are enclosed in dot- 
ted black parallelograms. Five branches are shown, with colours 
and line-styles indicated. The solid gray lines surrounding the 
branches enclose half the mass. Of the 10 potential branch pairs, 
only five pairs interact within <i cr jt = 200 kpc h . Of those, all 
result in a merger, except for the C-D pair. Branch E is tidally- 
strippcd below resolution before merging on to branch D (an open 
circle on branch E at z ~ 1.2 is missing). 



coordinates) of the haloes identified in the simulation. The 
radii are proportional to the logarithm of the mass. The solid 
gray lines flanking the branches enclose half of the mass of 
each halo, giving a rough sense of the 'size' of these ob- 
jects. The solid red circles denote mergers and the black 
dotted rectangles (or parallelograms) represent close inter- 
actions. The branches are said to be 'interacting' if they 
are within some critical distance d cr it of each other (in the 
Figure, d crit = 200 kpc /i _1 ). 

The merger history occurs as follows. First branch A 
(dark-green dot-dashed lines) merges with branch B (blue 
long-dashed lines). Of these two, only branch B continues 
on to later epochs after the merger. Another way to put 
this is that branch A 'merges on to' branch B because, right 
before the merger, the halo on branch B is more massive than 
the halo on branch A. Later on, branch C (magenta short- 
dashed lines) also merges on to branch B. This is followed by 
branch B merging on to branch D (dark-brown solid lines). 
Lastly, branch E (dotted long-dashed green lines) merges on 
to branch D. In this merger tree, branch D is the principal 
branch of the final halo rooted at z = 0. 

Not all branches begin at the same redshift. For in- 
stance: branches A and B appear simultaneously. Later, 
branch D appears, followed by branch C. Branch E is born 
long after most of the other branches have already merged. 
Physically, the haloes might have been around from much 
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Figure 2. The mass growth of two interacting haloes. The more 
massive of the two begins at ~ 1.2 X lO 12 M0h — 1 and grows 
up to ~ 5 X 1O 13 M0/i~ 1 (solid light-red line). The secondary 
halo begins at ~ 3 X 10 10 Mq}i~ , reaches a maximum mass of 
~ 3 X IO^MqIi- 1 (at z ~ 1), and finishes at z = with ~ 
2.5 X l0 10 M e h- 1 (dashed dark -green line). The horizontal blue 
dotted lines represents the threshold that selects which branches 
make it to the catalog (set at 100 particles = 8.6 X 10 10 Mq}i~ 1 ). 
Branches are selected by whether they cross this threshold at 
some point in their history, regardless of their final mass. 



Moreover, branch E appears after branches A, B and C have 
already merged on to branch D. In other words, branch E 
never 'coexists' with branches A, B and C. Although there 
are five closely-interacting pairs, only four end up as merg- 
ers. Pair C-D interacts briefly but eventually gets separated 
beyond d cr it- Branch C never has the opportunity to merge 
on to branch D, and instead merges on to branch B. It is 
often tempting to assume that all close pairs merge eventu- 
ally. One of the goals of this paper is to quantify the fraction 
of close pairs that never merge. 

Figure[T]only shows a single merger tree. However, many 
massive haloes often have plenty of substructures. The ap- 
propriate thing to do in such situations is to follow the 
merger histories of both the central dominant halo and the 
satellites. In this paper, I pair up all branches in a given 
group, even if they belong to different merger trees within 
that group. Often I will use the term 'never-merging' to refer 
to pairs that do not merge because either one of the members 
in the pair gets absorbed by a third halo or because both ha- 
los survive as separate substructures of a larger background 
halo until the present time. 

Lastly, one could also pair up branches belonging to 
different groups; i.e., pairs of haloes that never end up grav- 
itationally bound to each other. These events could be im- 
porta nt if one tries to u se close pairs as a proxy of mergers 
(e.g.. IWetzel et all 120081 ). especially if large separations (of 
the order of 1 Mpc hT 1 ) are used. This contribution is more 
likely to be at the few-percent level if smaller separations 
(of the order of 200 kpc h -1 ) are used (private communi- 
cation with the reviewer). For the sake of simplicity, in this 
work I ignore those pairs whose haloes end up in different 
host haloes. The effects of this particular limitation will be 
addressed in future work in the context of close-pair counts. 



earlier epochs. However, the simulation only identifies haloes 
with at least 20 particles. 

The growth of haloes within larger haloes (i.e., sub- 
haloes) is commonly thwarted by the rich environments 
surrounding them, whic h tidally str i p them of t heir dif- 
fuse outer layers (e.g.-. Moore et alj 1 19961. 19991; Tormen 



19971 : iTormen et al.lll998l: iGhigna et al.lll998l : lKivDin et al 
19991 : iTaffoni et al.ll2003l : lGiocoli et al.ll2008l ) In many cases 



haloes go below the resolution of the simulation, and are no 
longer identified. Such is the case of branch E, which disap- 
pears before merging on to branch D. Extra care must be 
taken when dealing with these situations. Namely, one must 
pay attention and make sure that both haloes are identified 
at every snapshot in the interaction. In this work I impose 
an extra condition on the mass which deletes many (but not 
all) of the small sub-haloes that get stripped apart by their 
hosts (see Section [2.4.11 below). With this condition, much 
less than one percent of the pairs have one member going 
below resolution (situations where both members go below 
resolution are non-existent under this condition) . For the re- 
mainder of the paper, these effects will be ignored because 
they happen so infrequently. 

A merger tree with N branches has N(N — l)/2 candi- 
date pairs. However, not all of them qualify as interacting. 
For instance, the merger tree in Figure [T] has 5 branches and 
5(5— 1)/2 = 10 candidate pairs; but 5 of those 10 pairs never 
count as 'close' pairs. Namely, the distance between branch 
A and either branches C or D is always greater than d C rit- 



2.3 The Stages of Evolution 

In this paper, the evolution of a pair of closely interacting 
haloes is described in terms of the following five stages: 

• Initial: When the two haloes both exist for the first 
time (above the resolution of the simulation), regardless of 
whether they are closely interacting or not. 

• Entry: When the two haloes are, for the first time, at a 
distance less than d cr it. This may occur either at the initial 
stage or later. 

• Closest: When the two haloes are both identified as 
separate entities by the simulation and are at their closest 
distance. 

• Final: The last time the two haloes are both identified 
by the simulation, regardless of whether they are closely 
interacting or not. 

• Fate: What ultimately happens to the pair: either (1) 
they merge, (2) one or both haloes are absorbed by a third 
halo or (3) both haloes survive until today. 

It is useful to clarify a few subtleties in these definitions 
with specific examples from FigureQ] First of all, notice that 
both haloes in the A-B pair appear simultaneously and are 
both very small (just above resolution, with mass ratio very 
close to one) at the initial time. In other cases, however, 
one of the two haloes may have some time to grow before 
its partner shows up in the simulation (pairs B-C and C-D, 
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and especially D-E - the latter with a mass ratio particularly 
different from unity). 

When the initial distance of the pair is below d C rit, the 
entry and initial times coincide (pairs A-B, C-D and D-E). 
For pairs A-C and A-D, on the other hand, the initial dis- 
tance is larger than d C rit; and the entry time occurs much 
later than the initial time. When the entry time does not 
correspond to one of the snapshots in the simulation (e.g., 
as in pairs B-C and B-D), a linear interpolation in log(l + z) 
is performed. 

The closest distance is often equal to the final distance, 
especially right before mergers (pairs B-C, B-D and D-E). 
In some cases, though, the closest distance does not match 
the final distance. Examples of these include pairs A-B and 
C-D; where the former results in a merger and the latter 
does not. 

The final epoch refers to the last time both haloes are 
identified in the simulation as separate entities. They need 
not be within d C rit: for instance the final time for pair C- 
D takes place at z ~ 2.7, with separation ~ 2d cr it- In this 
particular case, it is interesting to know when the pair is 
within ddit for the last time (i.e., z ~ 4 for pair C-D). In 
Section 13.51 I use this information to estimate the typical 
'lifetime' of interactions. 

The ultimate fate stage is formally not considered part 
of the evolution of the pair. Failing to make this distinction 
may lead to over-counting the number of instances a halo 
participates in an interaction. For example, the halo result- 
ing from the merger of branches A and B participates in both 
pairs B-C and B-D, but not in A-B itself. (Once the merger 
has happened, the single remnant halo cannot be counted 
as part of its own halo pair!) It is important to probe the 
interval between the final and the fate time. If this portion 
is not included towards the lifetime of the pair, one could 
conclude that pair D-E has zero duration, whereas in reality 
it lasts from z ~ 2.4 down to z ~ 1.8 (slightly less than 
three full snapshots). Nevertheless, as mentioned before, in 
this paper I impose a mass condition that deletes most pairs 
involving haloes going below resolution (see section 12.4.11 ) 
- thus situations like the D-E pair are extremely uncommon 
with this condition. Lastly, if the final epoch takes place to- 
day, the fate epoch is set equal to the final epoch. This is the 
case for pairs involving different members of a group today. 

2.4 Constructing the Catalogue 

In this Section I explain the construction of the halo pair 
catalogue used in this paper, which can be understood as a 
series of selection criteria. I defer discussions on the adopted 
choices and assumptions to Section f4. II 

2.4-1 Mass Selection 

First, I only consider branches that reach a critical mass 
M cr it at least once at some point throughout their life- 
time. Isolated haloes typically reach their maximum mass 
at 2 = 0. However, recall that a halo immersed in a larger 
group (i.e., a subhalo) may be substantially dissolved by the 
gravitational tidal field of its host. In this case, it is possible 
that at some point a given halo had mass > M cr it, but ends 
up with mass < M cr it. In some extreme cases, it may even go 



below the resolution of the simulation. It must be mentioned 
that for general branches (i.e., those not obeying the above 
mass condition), this situation are not uncommon. Never- 
theless, once the above mass condition is imposed, such sit- 
uations become very infrequent (less than one percent of 
the pairs have one member going below resolution at the 
late stages of interaction). 

In this work, massive branches are kept even if they 
lose a substantial fraction of their mass at later epochs. 
To select these branches, the appropriate thing to do is to 
choose them by their maximum mass, not their final mass. 
Figure [5] illustrates this situation by showing the mass his- 
tory of two interacting haloes that are both accepted into 
the catalogue by this criterion. The horizontal dotted blue 
line denotes the mass selection threshold (in the Figure, 
M cr u = 8.6 x IO^Mq/i -1 = 100 x TO P ). The solid light-red 
line represents the more massive halo, while the secondary 
halo is shown as a dark-green dashed line. The first branch 
grows from ~ 1.2 x IO^Mq/i" 1 up to ~ 5 x IO^MqTT 1 . For 
the most part, the mass of this halo increases monotonically 
(except for the brief dip near z ~ 0.2). At all times this 
branch is above the selection threshold. The second branch 
begins at ~ 3 x 10 10 ' Mq^ 1 , below the selection threshold. 
At z ~ 3.5, it crosses the threshold - reaching a maximum 
mass of ~ 3 x 10 11 Mq/i -1 at z ~ 1. At that point, it begins 
to be tidally stripped of its mass and goes below the selec- 
tion threshold at z ~ 0.7. By z = 0, it ends up with a final 
mass of ~ 2.5 x 10 11 M Q /i _1 . If the selection were based on 
the mass today, this branch would have been thrown out, 
despite the fact that at some point it had a mass four times 
greater than the threshold! 

In summary, only branches that reach a critical mass 
M cr it at some point in their existence are paired up in the 
catalogue. This maximum mass is usually attained at the 
end of their lives as isolated haloes, right before they are 
accreted by larger haloes. This is often called the 'infall 
mass', and many authors consider it an important ingre- 
dient in galaxy for mation models (e.g.. IConrov et aUl2006l ; 
iMoster et al]|201Ch . 

2.4-2 Co-existence 

With the list of branches satisfying the above mass selection 
criterion, the next step is to collect pairs. A tree with N such 
branches can have up to N(N—l)/2 pairs. However, in many 
cases, a branch disappears (it merges on to another branch) 
long before the candidate partner shows up in the simulation 
(see discussion about branch E in regards to branches A, B 
and C; in Section [2T2l and Figure 0. The correct way to check 
whether two branches 'co-exist' is to see if the intervals from 
birth (when the branch first appears in the simulation) to 
fate of the two branches ever match. Special care must be 
taken when branches go below resolution, typically towards 
the end of their existence. 



2-4-3 Distance Selection 

Once the above two criteria have been satisfied, one can 
readily identify the initial and final stages, as explained in 
Section 12.31 The relative distance can be calculated for all 
those snapshots for which both haloes are identified. With 
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Set 


<^crit 


Mcrit/MQh- 1 


Merit /m p 


Number 


C 


1 Mpc h' 1 


8.6 x 10 10 


100 


502,705 


vc 


200 kpc h" 1 


8.6 x 10 10 


100 


32,582 


MC 


1 Mpc h- 1 


8.6 x 10 11 


1,000 


13,637 


MVC 


200 kpc h' 1 


8.6 x 10 11 


1,000 


2,067 



Table 1. Halo pair sets shown in the panels of all Figures in 
Section [3] The abbreviations are: C (Close), VC (Very Close), 
MC (Massive Close) and MVC (Massive Very Close). Distances 
are in comoving coordinates. Recall that m p = 8.6 X 1CPMqH~ 1 . 
For reference, C and MC are the Close Sets, VC and MVC are 
the Very Close Sets, C and VC are the non-Massive Sets, and 
MC and MVC are the Massive Sets. 



few up to 10 Mpc ft -1 , but never beyond 30 Mpc ft . 
There is, however, a small fraction of pairs (~5%) with 
^initial < dcrit. The most common closest separation is near 
dcrit- Recall that close branches are selected so d c ioscst < 
d cr it, which is why this distribution is truncated at d cr it- 
At the final stage, a large fraction of pairs end up at dis- 
tances beyond d cr it- For the rest, the final-distance distribu- 
tion resembles the closest-distance distribution. In all cases, 
a non- negligible fraction of pairs never merge (labeled as 
' fate with two haloes'), and many of those even end up at 
distances greater than d cr it- Below, I describe these trends 
for each stage separately. 



this, the next step is to find the minimum distance dmm 
(the closest stage, Section T2.3p . and check if the dmm < dcrit 
criterion is met. 

244 The Halo-Pair Sets 

The following sets of halo-pairs are considered: the Close 
Set (C), the Very Close Set (VC), the Massive Close 
Set (MC) and the Massive Very Close Set (MVC). These 
sets are defined by the following critical parameters: 

• C: (dcrit, M crit ) = (1 Mpc ft -1 , 8.6 x 10 10 M qj h- 1 ). 

• VC: (d C rit, M C rit) = (200 kpc ft" 1 , 8.6 x lO^Moft" 1 ). 

• MC: (d C rit, M C rit) = (1 Mpc ft" 1 , 8.6 X 10 11 M Q h' r ). 

• MVC: (dcrit, Merit) = (200 Mpcft- 1 ,8.6xl0 11 M Q ft- 1 ). 

Notice that the latter three sets are simply sub-catalogues 
of set C. Often in the text, and in all the remaining Fig- 
ures, I use the following terminology: C and MC are the 
Close Sets (with d cr it = 1 Mpc ft" 1 , upper panels); VC and 
MVC are the Very Close Sets (with d cr it = 200 kpc ft -1 , 
lower panels); C and VC are the non-Massive Sets (with 
Merit = 8.6 x 10 10 M©ft- 1 , left-hand panels); and MC and 
MVC are the Massive Sets (with A/ C r it = 8.6x lO n M ft -1 , 
right-hand panels). For quick reference, Table[T]indicates the 
assigned values of d cr it and M cr it, and the above terminology. 



3 RESULTS 

In this Section I present the distribution of several quantities 
as a function of what stage in their evolution the pairs hap- 
pen to be. I also address relevant lifetimes and the effects of 
modifying d cr it and M cr it. Only results are reported here. I 
postpone associated interpretations and discussions for Sec- 
tion(4] I do this for the sake of clarity because the forthcom- 
ing discussion requires having described all the Figures in 
this Section already. 

3.1 Relative Distances 

Each panel in Figure [3] shows each of the four sets in Table [T] 
The vertical line is drawn at d cr it to indicate the separation 
at the entry stage. Relevant colours and line-styles are indi- 
cated in all Figures. 

Typically, pairs have large initial distances, between a 



3.1.1 Initial Distances 

The initial-distance distribution for the Close Sets peaks at 
~ 10 Mpc ft -1 , is narrow, and drops sharply at low sep- 
arations. For the Very Close Sets, the distribution peaks 
around 8 and 3 Mpc ft -1 (for the VC and MVC sets respec- 
tively), is wider, and has a small secondary peak at ~ 100 
kpc ft -1 before dropping sharply at low separations (the 
secondary peak is more evident for the MVC set). Reducing 
dcrit does not alter the distributions significantly, except for 
the appearance of a bump slightly below d cr it and a trough 
slightly above this value. Increasing M cr it also enhances the 
low-separation regime, but to a very small degree. I must 
emphasize that the initial stage depends strongly on the 
resolution of the simulation. Nevertheless, it is interesting 
to show the highest-redshift status of the close pairs one 
can study with the Millennium Simulation. 



3.1.2 Closest Distances 

For the most part, the closest-distance distribution increases 
monotonically with increasing distance and is truncated at 
dcrit- For the Very Close Sets, the peak coincides with the 
location of the secondary peak in the initial-distance dis- 
tribution. Reducing d cr it strongly increases the contribution 
from pairs with small minimum separations (in particular, 
below 100 kpc ft -1 ). Increasing M cr it leaves the distribu- 
tion almost intact, except for a slight enhancement at small 
separations. 



3.1.3 Final Distances 

The shape of final-distance distribution mimics that of the 
closest-distance distributions up to d cr it, and then decreases 
beyond that (with an extra small bump at ~ 800 kpc ft -1 
for the Very Close Sets) . The final distance never exceedes 5 
Mpc ft - . For the Close Sets, the final distance distribution 
peaks at ~ 800 kpc ft -1 , while the peak is at ~ 200 kpc ft -1 
for the Very Close Sets. 

Decreasing d cr it shifts the distribution to smaller sepa- 
rations, enhancing the number of small-distance pairs (with 
dfinai < 200 kpc ft -1 ) quite significantly. However, this 
change is not so strong for pairs with dfl na i < d cr it- On the 
other hand, the contribution to other regimes (> 200 kpc 
ft -1 ) is suppressed with decreasing d cr it- The effects of in- 
creasing Afcrit are similar to those of decreasing d cr it, but 
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Figure 3. The relative-distance distribution of halo pairs. The set of halo pairs presented (see Table [T} is indicated in each panel: C 
(upper-left), VC (lower-left), MC (upper-right) and MVC (lower-right). The colours and line-styles are indicated in the Figure. The vertical 
line represents ci cr jt. Only pairs that never merger are presented for the fate stage (indicated by the 'two haloes' label). For reference, 
the text uses the following definitions to describe the sets: Close Sets (d cr it = 1 Mpc h~ 1 ; upper panels), Very Close Sets (c( cr i t = 200 
kpc /i -1 ; lower panels), non-Massive Sets (M cri t = 8.6 X 10 10 Moft, _1 ; left-hand panels), and Massive Sets (M cri t = 8.6 X 10 11 M ft, _1 ; 
right-hand panels). See Section 12.31 and 12. 4. 41 respectively for definitions and descriptions of the stages and sets. 



milder. However, while decreasing d cr i t enhances the sec- 
ondary bump present in the Very Close Sets, increasing M cr it 
actually suppresses the bump, but only slightly. 

3.1.4 Fate Distances (Two Haloes) 

I only present the fraction of fate distances of pairs that 
never merge. (Recall that never-merging pairs are those that 



either get split by a third party or survive as two entities 
until today.) For the fate two-halo distribution, I will of- 
ten speak of the 'absolute' and 'relative' contribution from 
a given regime. The former is normalised to the total pop- 
ulation of pairs, while the latter is normalised to the subset 
of pairs that never merge. In other words, the sum of all 
relative contributions is 100%, while the sum of all absolute 
contributions is one of the four percentages quoted above. 
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Figure 4. The redshift distribution of halo pairs. The set of halo pairs presented (see Table [TJ is indicated in each panel: C (upper-left), 
VC (lower- left), MC (upper-right) and MVC (lower-right). The colours and line-styles are indicated in the Figure. Only pairs that never 
merger are presented for the fate stage (indicated by the 'two haloes' label). For reference, the text uses the following definitions to 
describe the sets: Close Sets (d cr ;t = 1 Mpc h~ 1 ; upper panels), Very Close Sets (e£ cr it = 200 kpc h —1 ; lower panels), non-Massive Sets 
(M crit = 8.6 X lO lo M0/i _1 ; left-hand panels), and Massive Sets (Merit = 8.6 X 10 n Mgh -1 ; right-hand panels). See Section l2~3l and 
12.4.41 respectively for definitions and descriptions of the stages and sets. 



The shape of the fate-distance distribution is similar to 
that of the final-distance distributions for the Close Sets, and 
to the secondary bump in the final-distance distributions for 
the Very Close Sets. Reducing d cr i t decreases the absolute 
number of pairs that never merge. It shifts the distribution 
to lower distances (the peak moves from ~ 800 — >-~ 400 kpc 
h^ 1 for the non-Massive Sets, and from ~ 800 — »~ 300 kpc 
h~ x for the Massive Sets). However, the relative fraction of 



pairs with df ate > d C rit increases. In other words, there is a 
shift to smaller separations; but this shift is actually to larger 
separations in units of d C rit- Increasing M cr it also suppresses 
the absolute contribution of pairs that never merge. This 
effect is more evident than that produced by reducing d cr it- 
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3.2 Redshifts 

Figure [4] shows the distributions of halo-pair redshifts at the 
various stages of their evolution. The initial- redshift distri- 
bution is dominated by pairs at very early epochs (z > 4), 
with some low-redshift contributions for the Massive Sets. 
The subsequent stages of interaction (i.e., entry —5- closest 
— > final ) happen at much lower redshifts (compared to the 
initial distribution), and are increasingly less frequent with 
high redshift. As the interaction progresses, the distribu- 
tion shifts to lower redshifts (the final distribution exhibits 
a small upturn at very low redshifts). This trend continues 
at the fate stage for those pairs that never merge. More than 
~50% of the never-merging fate pairs occur at z < 0.5, while 
less than ~10% take place at z > 1. 

Reducing d cr i t only enhances the fraction of initial pairs 
below z ~ 1 slightly. This change makes the entry-redshift 
distribution shallower in the z < 3 regime, and introduces a 
sharp drop afterwards, bringing its shape more in line with 
the closest and final-redshift distributions. The same flatten- 
ing and drop appear for the closest and final distributions. 
The fate distribution is suppressed more strongly as a func- 
tion of redshift, increasing the discrepancy with the other 
three distributions described here. Increasing M„it also en- 
hances the fraction of initial pairs below z = 1. Additionally 
the peak near z ~ 6 becomes more narrow, and a trough at 
z ~ 3 appears. For the entry distribution, the high-redshift 
regime is slightly suppressed (but not completely erased) for 
the Close Sets, and completely cut off beyond z > 4 for the 
Very Close Sets. The closest and final distributions are also 
cut-off completely beyond z > 4. The fate distributions are 
completely suppressed beyond z > 3, and enhanced in the 
z < 0.5 regime. 



3.3.1 Initial Total Masses 

Reducing d cr it does not alter the initial-mass distribution 
by much, except for a small enhancement in the large-mass 
regime (with masses > a few x10 13 Mq/i _1 ). A similar ef- 
fect happens when M cr i t is increased. This distribution is 
strongly determined by the resolution of the simulation, and 
its features will not be addressed in much detail. Neverthe- 
less, it serves as reference to see how much the branches 
grow during the interacting stages. 

3.3.2 Entry, Closest, and Final Total Masses 

In this Section I describe the entry, closest and final stages 
together because they follow similar total-mass distribu- 
tions. These distributions are broader than the initial-mass 
distribution, and peak at a higher values: slightly above 10 11 
and IO^Mq/i" 1 for the non-Massive and Massive Sets re- 
spectively. Relative to the initial stage, the typical total mass 
increases by a factor of ~ 2 for the non-Massive Sets, and 
by a factor of ~ 20 for the Massive Sets. Also, the large- 
mass contribution (with masses > a few xlO 13 MqIi^ 1 ) is 
no longer negligible. 

Reducing d cr it does not alter the distributions much, ex- 
cept that the high-mass drop becomes shallower, enhancing 
the contribution from that regime. This is particularly no- 
ticeable for the Massive Sets. Increasing M cr i t shifts the peak 
of the distributions from ~ 10 11 -> 10 12 M Q ft _1 . This sup- 
presses the small-mass regime, and enhances the large-mass 
regime. The location of the new peak coincides with that of 
the small bump found in the initial total-mass distribution. 
The large-mass plateau is enhanced more significantly here 
than when d cr i t is reduced. 



3.3 Total Masses 

Figure [5] shows the distribution of the sum of the two masses 
(hereafter, the total mass) in a halo-pair at various stages 
of their evolution. At the initial stage, the majority of the 
total masses are rather small (peaking below 10 ll Mgft — l , 
with a negligible contribution from masses greater than a 
few x 1O 13 M0/i _1 ). The transition from the initial to the 
entry stage populates the large-mass regime, evacuating the 
small-mass regime as a result. This shift towards more mas- 
sive pairs is more evident for the Massive Sets. As the in- 
teraction progresses (i.e., entry, closest and final stages), 
the very large-mass regime (> a few xlO 13 Mq/i" 1 ) con- 
tinues to grow, while the rest shrinks. (This transition mass 
scale is not very evident for the MVC set.) These shifts are 
very small relative to the previous initial — > entry shift, but 
not entirely insignificant. The fate -mass distribution of pairs 
that never merge does not resemble the others because it is 
controlled by third parties that absorb the pairs, not by 
the pairs themselves. Notice that these pairs are incorpo- 
rated preferentially into extremely massive haloes (ending 
up with masses typically greater than 10 14 Mq/i _1 ). These 
situations most likely involve two subhaloes interacting for 
a while until they are split apart by the central halo. 



3.3.3 Fate Total Masses (Two Haloes) 

Figure [5] shows the fate total-mass distribution, which is 
completely different from the distributions at the other 
stages. This distribution increases monotonically with mass, 
except at the very high masses, where a small dip is present. 
Reducing d cr it and M CI - lt both have the same effect: to 
enhance the small-mass end and suppress the large-mass 
regime. This effect is stronger for decreasing Merit than for 
decreasing d C rit- Notice that for the non-Massive Sets, the 
peak of the fate-mass distribution is three orders of magni- 
tude higher than that of the final-mass distribution, while 
only two orders of magnitude for the Massive Sets. 

3.4 Mass Ratios 

Figure [6] shows the distribution of the mass ratios of the 
haloes in the pairs at the various stages of their evolution. In 
this work, the mass ratio is defined as the ratio of the small 
to the large mass in a given pair (i.e., this quantity is always 
^ 1; small values correspond to disparate pairs while large 
values correspond to nearly-equal pairs). In particular, the 
primary versus secondary ranking of the haloes in the pairs 
need not remain fixed at all times: at times the haloes may 
switch roles. For instance, a halo may begin more massive 
than its companion, but then the companion could grow 
quickly for some period and quickly surpass its partner. In 
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Figure 5. The total mass distribution of halo pairs. The set of halo pairs presented (see Table[TJ is indicated in each panel: C (upper-left), 
VC (lower- left), MC (upper-right) and MVC (lower-right). The colours and line-styles are indicated in the Figure. Only pairs that never 
merger are presented for the fate stage (indicated by the 'two haloes' label). For reference, the text uses the following definitions to 
describe the sets: Close Sets (d CT n = 1 Mpc h~ 1 ; upper panels), Very Close Sets (rf cr jt = 200 kpc h —1 ; lower panels), non-Massive Sets 



8.6 X 10 lu M Q /i _1 ; left-hand panels), and Massive Sets (M c 



12.4.41 respectively for definitions and descriptions of the stages and sets. 



8.6 X 1O 11 M /i _:l ; right-hand panels). See Section l2~3l and 



these situations, the ranking of the mass ratios are adjusted 
accordingly. 

The resolution of the simulation, and the criterion that 
branches must surpass a critical mass threshold, may affect 
the mass ratio distribution (especially for the initial stage, 
see Section |4] for a discussion). Nevertheless, these results 
are worth presenting because (1) they provide additional 
clues describing how the pairs in the simulation grow, and 



(2) they help us understand the effects of changing d cr i t and 
M cr it more clearly. 

The initial , entry, closest and final mass-ratio distribu- 
tions all increase monotonically with mass ratio. The distri- 
bution is sharply peaked initially at very high mass ratios 
and becomes broader as the pairs evolve (initial — > entry 
— ► closest —> final). During this process, the low mass-ratio 
regime gets gradually more populated. In other words, the 
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Figure 6. The mass-ratio distribution of halo pairs. The set of halo pairs presented (see Table[T} is indicated in each panel: C (upper- 
left), VC (lower-left), MC (upper-right) and MVC (lower-right). The colours and line-styles are indicated in the Figure. Only pairs that 
never merger are presented for the fate stage (indicated by the 'two haloes' label). For reference, the text uses the following definitions 
to describe the sets: Close Sets (d cr ; t = 1 Mpc h —1 ; upper panels), Very Close Sets (d cr it = 200 kpc lower panels), non-Massive 

Sets (Merit = 8.6 X 10 10 Afo/i _1 ; left-hand panels), and Massive Sets (M crit = 8.6 X 1O 11 M /i _1 ; right-hand panels). See Section EOl 
and 12.4.41 respectively for definitions and descriptions of the stages and sets. 



masses in pairs tend to become more discrepant with time. 
Notice that as this evolution progresses, the distribution is 
suppressed at mass ratios > 0.1 and enhanced at mass ra- 
tios < 0.1. Lastly, the closest and final mass ratios have 
nearly identical distributions, while the entry distribution 
is similar to these two, except for a small difference in the 
low mass-ratio regime. As in the analysis of the total mass 
(Section 13. 3|) . the fate mass-ratio distribution is completely 



different from the other stages. For instance, the fate mass 
ratios are preferentially < 0.1 (recall too that the sum of the 
fate-halo masses is preferentially large). 



3.5 Lifetimes 

In this Section I estimate the typical duration of close inter- 
actions. The following relevant time intervals are considered: 
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Figure 7. The interaction-duration distribution of halo pairs. The set of halo pairs presented (see Table [TJ is indicated in each panel: 
C (upper- left), VC (lower- left), MC (upper-right) and MVC (lower-right). The colours and line-styles are indicated in the Figure. The 
total duration refers to the time from the initial until the fate stage. The approach time goes from the initial until the entry stage. The 
interaction lasts from the first until the last time the pair is identified within d cr it- For reference, the text uses the following definitions 
to describe the sets: Close Sets (d C rit = 1 Mpc ft -1 ; upper panels), Very Close Sets (c£ C rit = 200 kpc lower panels), non-Massive 

Sets (M CT it = 8.6 X 10 10 A/ O h- 1 ; left-hand panels), and Massive Sets (M crit = 8.6 X 10 11 Af Q /i- 1 ; right-hand panels). 



• Total Time: The interval from the initial to the fate 
time. Here, unlike in previous discussionsSection, fate here 
is not limited to the non- merging pairs. 

• Approach Time: The interval from the initial to the 
entry time. This refers to the approach time from large dis- 
tances down to the first time the pair is within d cr it- Notice 
that those pairs with (initial < d C rit have zero approach time. 



• Interaction Time: The interval from the first to the 
last time the pair is within dcrit. If dmitial < dcrit, the inter- 
action begins at the initial time. If dtnitial > dcrit, it begins 
at the entry time. If df ate < d cr it (i.e., mergers or sub-critical 
non- merged pairs), the interaction ends at the fate time. If 
dfate > dcrit, the last time at which the relative separation 
crosses the d cr it threshold is sought. If this time does not co- 
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incide with one of the snapshots in the simulation, a linear 
interpolation in log(l + z) is performed. 

It is important to mention that the approach time de- 
pends on the initial time, which itself is strongly determined 
by the resolution of the simulation. This issue also affects the 
estimate of the total time. Despite this limitation, it is in- 
teresting to study the typical timescales experience by close 
pairs in the Millennium Simulation, and the effects of mod- 
ifying d cr it and Merit- The interaction time is less suscep- 
tible to these resolution effects because the mass-criterion 
adopted in this paper deletes most of the branches that are 
most prone to go below the mass resolution of the simulation 
more quickly due to tidal stripping. Nevertheless, additional 
improvements using dynamical-friction recipes can always 
be incorporated in future works (which will most likely in- 
crease the duration of the interaction times). The results 
presented here serve as a reference for such developments 
(see Section [4] for additional discussions on these issues). 

The three lifetime distributions are shown in Figure 
The total -time distribution peaks sharply at long durations, 
with over 80% with ttotal > 5 Gyr (and over 90% with such 
long durations for the Massive Sets). The approach-time 
distribution resembles the total-time distribution, but with 
some contributions from shorter durations. Over 65% of the 
pairs have t apP roach > 5 Gyr (and over 80% have such long 
durations for the Massive Sets). The interaction -time dis- 
tribution is dominated by durations between approximately 
0.2 — 4 Gyr. The distribution has a long tail at low-durations 
(with f interaction < 0.1 Gyr), especially for the Very Close 
Sets (~20% for the VC set, and ~12% for the MVC). 

Reducing either d cr it or M cr it does not affect the distri- 
butions significantly. There is a small enhancement in the 
intermediate time regime (approximately 1 — 8 Gyr), but 
not for the longest-time regimes. The effect of changing d C rit 
is slightly stronger than that caused by changing Merit- In- 
creasing Merit shifts the lower-end of both the total and the 
approach-time distributions to higher values. The most ob- 
vious change is that adopting a smaller d cr i t makes short 
interactions more common, while choosing a larger M cr it 
makes the approach-times longer. Lastly, when both modifi- 
cations are applied, interactions lasting longer than ~8 Gyr 
are essentially absent. 



4 DISCUSSION 

This Section discusses the assumptions and provides inter- 
pretations for results presented in this paper. A description 
of related work is included. 



4.1 The Assumptions in the Catalogue 

Before analysing and intepretting the results of Section 
I devote this section to discuss the assumptions adopted in 
this paper. 

The catalogue analysed in this work consists of halo 
pairs that (1) come within some comoving distance dcrit and, 
(2) in which both haloes reach a mass greater than some 
threshold value M cr it at some point in their history. 

The primary (largest) catalogue is dubbed the 'Close 
Set', and it has d C rit = 1 Mpc h -1 and M cr it = 8.6 x 



10 1 Mq/i -1 (equivalent to 100 simulation particles). Sub- 
catalogues are then selected from this set by either decreas- 
ing the critical distance to d cr i t = 200 kpc h~ 1 or / and 
by increasing the critical mass to Merit = 8.6 x 10 1 Mq/i~ 
(equivalent to 1,000 simulation particles). From these com- 
binations, the following additional sets arise: the Very Close 
Set, the Massive Close Set and the Massive Very Close Set. 
See Table [T] for a summary. One of the goals of this work is 
to quantify the effects of modifying these selection criteria. 

The larger critical distance adopted for the primary cat- 
alogue (dcrit = 1 Mpc h^ 1 ) is not as large as it appears be- 
cause it is in comoving coordinates. In these coordinates, this 
threshold is compara ble to the large- s epara tion tail of bi- 
nary quasars found by I Foreman et. al] (2009). Indeed, these 
authors report 38 binaries (out of 85) with comoving sepa- 
rations larger than 200 kpc h^ 1 , and 3 binaries with separa- 
tions larger than 1 Mpc h~ \ One could argue that such pairs 
are not truly associated. However, these authors found their 
measurements to be enhanced relative to simple extrapola- 
tions of the quasar correlation function, signaling that they 
may indeed be physically influencing each other. 

The distance criterion used here is still larger than some 
of the popular thresholds found in the literature (< 50 physi- 
cal kpc are commonly used to select pairs with enhanced star 
formation). Nevertheless, keeping this exaggerated choice is 
valuable for a variety of reasons. First, it facilitates future 
studies of halo triplets and N-tuplets. Such investigations 
can shed light on the evolut ion of close compact groups 
- e.g ., like Stephan's quintet JStephanl 18771; Renaud et al.l 
2010) - and multip l e-quasar system s (Diorg ovski et al .1120071 ; 
Alonso et all 120081 ; iLiu et all [ioilal lbllcl). Similarly, the ef- 



fects of the local environment surrounding the pai r s can 
be tracked (e.g., iBarton et all 120071 ; lEllison et all |2010| . 
and references therein). Lastly, this large proximity choice 
will be useful for future studies of sub-Mp c clustering 

iMvers et all |200t1 ; IShen et all 



studi e s of quasars (e.g 
|2010| ; iBonoli et all 120101 ) 



and /or luminous red g a laxies 



200S 



(LRGs) (e.g., iBell et alj|2006l: iMasiedi et all I2OO 
Wake et al.1 l200fj, l200a Iwhite et all 120071; IConrov et al 



2007| ; lAlmeida et all |200S|; iBrown et al.1 120071: ICool et al 



20081 ; Ivan Dokkum et al.ll2010l ; iToieiro et al.ll201lh 

Alternatively, one could select pairs by whether they 
come within a distance below a fixed multiple of the sum 
of th e two virial radii (e.g., ISinha &i Hollev-Bockelmannl 
2011). However, the appropriate choice of separation for 
interaction- induced phenomena (e.g., star-f ormation en- 
hanc ement) may depend on environment (e.g.. lEllison et all 
2008). For the sake of simplicity, a selection based on centre- 
to-centre distance is preferred at this point. 

The mass criterion adopted in the primary catalogue 
(Merit = 100 x m p ) is introduced mainly for computa- 
tional convenience. Including all possible branches without 
any mass constraint is extremely cpu-intensive. Obviously, 
the problem of interacting dwarf galaxies is very interest- 
ing. However, using the Millennium Simulation for this pur- 
pose is not entirely appropriate, since modeling such galax- 
ies with less than a hundred particles is a very crude ap- 
proximation. For these systems, using something like the 
Millennium-II Simulatio n (with Lb ox = 100 Mpc h -1 and 
m p = 6.9 x 10 6 M m /i~ 1 . lBovlan-Kolchin et al-lkoOSl ) is more 
suitable. I reserve an extension of this catalogue in that di- 
rection for a future project. 
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4.2 Interpretation of the Results 

The analysis of Section [3] provides a relatively simple pic- 
ture that incorporates five basic stages: initial — >• entry — > 
closest — s> final — > fate. Below I use the wealth of statistical 
information gathered in Section [3] to make more sense of the 
details behind this process. 



4-2.1 The Initial Stage 

When both haloes are first identified (above the resolution 
of the simulation), they are typically found at distances of 
a few up to 10 Mpc ft -1 , but never beyond 30 Mpc h" 1 
(Figure [3} . Figure [3] shows that the pairs usually begin at 
redshift ~ 5 — 7, with only about 5% at z < 1 (except 
for the MVC set, which has a ^20% contribution from this 
redshift regime - see the foregoing discussion for a plausible 
interpretation). The majority of the total masses in the pairs 
are initially small: about 60% are in the (5-10)xl0 10 M Q /i- 1 
range, ~35% are in the (10 11 — 10 12 ) x MqH' 1 range, and 
only ~5% have masses greater than 10 12 MqH^ 1 (Figure^). 
About 90% of the initial mass ratios are greater than 0.1 
(Figure • 

The shape of the initial-distance distribution (Figure [3} 
depends mainly on two competing factors. As the size of a re- 
gion around a halo increases, more and more neighbours are 
within reach. Eventually, this quantity drops precipitously 
because the number of branches that participate in the as- 
sembly of the final host halo is finite. One would expect that 
decreasing the adopted d cr it would preferentially select pairs 
with smaller dinitiai- Overall, this is not the case, except at 
small distances: the Very Close Sets (lower panels of Fig- 
ure [3]) exhibit as a small secondary bump at 100 kpc 
and a trough at 300 kpc ft -1 . See below for a possible inter- 
pretation of this feature. Indeed, it is reasonable to speculate 
that with much smaller values of d C rit, this bump would be- 
come more prominent relative to the large-separation peak. 

A large fraction of initial halo pairs are 'major' (with 
mass ratio > 1/10, Figure This is no surprise. At 
early times, simulated haloes are expected to have masses 
slightly higher than the resolution limit (20 x m p = 1.72 x 
10 9 Mq/i -1 ), making the corresponding mass ratio of any 
two haloes < 1. In reality, one of the haloes often has to 
'wait' for the other to be born (e.g., branches C and D in Fig- 
ure [TJ. During this wait period, the primary halo has some 
time to grow. However, this growth cannot be too strong, 
nor too frequent - otherwise the contribution from pairs with 
low mass ratios would be more significant. 

The nature of the initial stage depends strongly on the 
capabilities of the simulation. A simulation with a smaller 
box would not generate massive enough pairs that begin at 
such early times, with such large separations. Likewise, it 
would not support enough massive final host haloes capa- 
ble of hosting such massive branches in large numbers. A 
simulation with higher resolution, on the other hand, would 
in principle shift the initial mass distribution to even lower 
masses and higher redshifts. However, since the majority 
of pairs begin at levels not too different from the resolu- 
tion threshold (which therefore consist, for the most part, 
of very similar haloes), the initial mass-ratio distribution is 
not expected to be very susceptible to this. 



4-2.2 The Entry, Closest and Final Stages 

Recall that the great majority of pairs fall from rather large 
distances. Only ~5% begin with dmitiai < dcrit- For the rest, 
it takes a few Gyr to go from dmitiai — ► dcrit- The entry, clos- 
est and final separations take place predominantly at low 
redshifts (~70-90% of these events happen at z < 1, while 
only ~0.5-5% occur at z > 3, Figure [4]). The total masses 
at these stages are a few xlO 11 Mo/i -1 for the non-Massive 
Sets, and a few xl0 12 M o fo _1 for the Massive Sets (Fig- 
ure As expected in hierarchical scenarios of halo growth, 
increasing M cr it shifts the masses to higher values, and red- 
shifts to lower values (left versus right-hand panels of Fig- 
ures [4] and [5] respectively) . 

Also, as in the initial stage, the mass-ratio distribution 
at these subsequent stages is dominated by large mass-ratios, 
but to a smaller degree (Figure [SJ. This is because by the 
time two haloes come within a critical distance, one of them 
might have had time to grow (e.g., pair D-E in Figure [TJ. In 
principle if a simulation with a much higher mass resolution 
were used - these distributions would include a less promi- 
nent contribution from low-mass ratios simply because there 
would be more pairs involving two small haloes at all red- 
shifts. Situations corresponding to massive haloes devouring 
minuscule ones would not be as important. This point high- 
lights the fact that mass-ratio statistics depends strongly 
on the resolution of the simulation. This limitation is not 
too problematic for the types of problems I intend to at- 
tack with this catalogue (e.g., interacting galaxies, binary 
quasars, etc.), where cuts in the type of haloes hosting gas- 
rich galaxies and supermassive black holes must be imposed. 

As the evolution of the pair progresses (entry — » clos- 
est — s> final), the corresponding redshift distributions shift 
towards more recent epochs. Also, the peak of the total 
mass distribution shifts slightly towards lower values. This 
behaviour is intriguing because one would expect mass to 
increase as the pairs (and the Universe) evolve. However, 
this shift toward lower masses takes place in the small-mass 
regime. Above a few xlO 13 M0/i -1 this trend is reversed: 
pairs get more massive as their evolution advances. (Al- 
though this effect is minimal for the MVC set: only a small 
fraction of low-mass pairs dissolve, and the rest remain in 
place for the duration of the three stages - see lower-right 
panel of Figure [5]) This evidence indicates that while mas- 
sive pairs follow their expected hierarchical-growth, smaller 
pairs are more susceptible to being stripped off of their mass, 
since most of them inhabit massive environments like clus- 
ters. A detailed analysis splitting this catalogue as a function 
of environment will shed light on this effect. This will be the 
subject of a future paper. 

Additional insight may be gained by looking at the be- 
haviour of the mass-ratio distribution as the evolution con- 
tinues. Recall that initially, ~90% of pairs are identified in 
the simulation with mass ratios > 1/10. By the entry time, 
this fraction drops to ~80% for the C set, ~ 70% for the MC 
and VC sets, and ~ 40% for the MVC set. The low mass- 
ratio regime, in turn, gets populated. In general, the masses 
of the haloes in pairs become more discrepant with time. In 
all cases this can be because, on average, either the primary 
haloes grows faster than the second one or the secondary 
halo dissolves more quickly. This effect is not expected to 
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be strong because, due to the resolution of the simulation, 
most pairs begin 'major' and remain more or less so. 

From the point of view of redshift, total mass, and mass 
ratio; the entry, closest and final stages behave rather simi- 
larly. The same cannot be said about the distance distribu- 
tion, making it a very powerful tool to learn about the dif- 
ferences of these stages. As in the initial case, distributions 
are expected to increase with distance (as more neighbours 
are seen) and drop (because groups have a finite number of 
branches). For the most part, the closest-distance distribu- 
tion only sees the rising portion. The exception is the MVC 
set (lower-right panel of Figure [3J, where the downturn is 
probed because groups have fewer massive branches satisfy- 
ing the smaller critical-distance condition. 

The distribution of final distances (along with the dis- 
tribution of never-merging fate distances discussed below) 
has a couple of interesting peculiarities. For sub-critical dis- 
tances, this distribution approximately follows the shape of 
the closest-distance distribution. These are the pairs that 
may, in principle, have the largest chance of merging. How- 
ever, there is also a large fraction of pairs ending up with 
dflnai > dcrit- This is expected because pairs are selected by 
whether or not they came within some critical separation, 
instead of by following imminent mergers back in time. 

Surprisingly, the fraction of super-critical final distances 
increases when d cr it is reduced (from ~45% —¥ ~57% for the 
non-Massive Sets and from ~39% — > ~45% for the Massive 
Sets - see portions to the right of the vertical line in Fig- 
ure[3]). A possible interpretation is as follows: as the adopted 
critical distance increases, the number of close neighbours 
a given halo can have also increases. While the number 
of massive neighbours increases, the number of less mas- 
sive haloes increases more substantially. In other words, the 
importance of having massive haloes in the vicinity dimin- 
ishes (relative to the less-massive neighbours) with increas- 
ing dcrit- Conversely, having a small d cr it means that the 
effects of third-party massive haloes in the neighbourhood 
(probably the main halo within the same group) become 
more dominant. One such effect is to split halo pairs: while 
two haloes may accompany each other for extended periods, 
one of them may be absorbed by a third massive halo be- 
fore they have the opportunity to merge. Another effect is 
to momentarily separate the pairs into very elongated orbits 
at the very last moments. Alternatively, these orbits might 
be elongated simply because when the distance between two 
haloes is comparable t o their virial radii, t heir orbits tend 
to be rather eccentric l|Zentner et al.ll2005i ). In either case, 
at the final stage, a large fraction of halo pairs are observed 
with separations much larger than dcrit- Indeed, it would be 
quite interesting to be able to track these last stages with 
more detail, using a simulation with much higher resolution. 

On a side note, if the importance of massive neighbours 
increases with decreasing d C rit, a larger contribution from 
low mass-ratio pairs should be expected . Figure [6] confirms 
this enhancement for all the stages mentioned here in (up- 
per versus lower panels). Nevertheless, this effect does not 
change the fact that the great majority of mass ratios remain 
more or less 'major' at all times. 

Finally, I must emphasize that the notions of closest and 
final stages should not be taken too literally. If two haloes 
merge, the closest and final distances they attain are strictly 
zero. My motivation to introduce these stages is not in this 



strict sense. Instead, my intention is to give an approximate 
picture (albeit snapshot-limited) of how the pairs evolve in 
the simulation. 



4.2.3 The Fate Stage 

It is important to keep track of the ultimate fate of the 
pair, even though this stage is not formally part of the evo- 
lution (although an infinitesimal of time before it actually 
is). As expected, reducing d cr it promotes the contribution 
of mergers, consequently reducing the fraction of never- 
merging pairs (from ~53% — ¥ ~42% for the non-Massive 
Sets, and from ~39% ->■ ~27% for the Massive Sets). 

For non-merging pairs, the fate-redshift distribution is 
strongly shifted towards recent epochs (fate-pairs with 2 < 
0.5 amount to ~50% for the non-Massive Sets, and ~60% 
for the Massive Sets - see Figure [4}. At these late epochs, 
massive groups and clusters become more prominent. These 
dense environments may be responsible for preventing the 
mergers of otherwise close pairs of haloes. 

Disruption of pairs can also be mitigated by adopting 
a larger Merit- This is because pairs of massive haloes are 
less prone to be disrupted by third-pary haloes. In principle, 
only super-massive neighbours could actually achieve these 
splittings, but such extreme objects are not very common. 

So far I have suggested that some pairs of haloes never 
merge because one or both haloes participating in the in- 
teraction are absorbed by a third party. In practice, none 
of the fate pairs get absorbed simultaneously by two dis- 
tinct haloes. Instead, it is common to have only one of the 
two haloes in the pair be eaten by a third object, leaving 
the other one wandering alone. (Moreover, in many cases 
the left-over halo is also absorbed by the same third ob- 
ject, most likely the central halo, but at a much later time.) 
This picture is supported by the following pieces of evidence. 
First, ~80% of the non-merging haloes end up with fate to- 
tal masses > 10 13 Mq/i _1 . Second, ~90% of this subset have 
fate mass-ratios < 1/100, supporting the interpretation that 
single, very massive third-party haloes are responsible for 
ripping the pairs apart. 

The fate distance distribution, just like the final dis- 
tance distribution discussed before, contains valuable spatial 
information. First of all, it is useful to compare the fraction 
of fate non-merging pairs with the fraction of super-critical 
final pairs. For the C, VC, MC, and MVC sets, recall that the 
super-critical final fractions are ~45%, 55%, 40% and 45%; 
while, on the other hand, the non-merging fate fractions 
(quoted above) are ~53, 42%, 39% and 27%. This demon- 
strates that the set of super-critical final pairs do not fully 
determine the set of never-merging fate pairs (read below 
for a possible interpretation). 

Nevertheless, if we focus on the super-critical fate 
regime instead, the corresponding fractions (relative to the 
general halo-pair population) are ~22%, 35%, 14% and 27%. 
This demonstrates two things. First, that the fraction of 
super-critical fate pairs is enhanced with reduced d C rit, dis- 
playing the same behaviour as their super-critical final coun- 
terparts. Second, that only a small fraction of the super- 
critical final pairs never merge. A significant fraction instead 
do return to sub-critical separations and merge eventually. 
This restriction of the snapshot-limited final stage can be 
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ameliorated by including velocity information across the last 
few snapshots. 

Likewise, the selection of close pairs by their closest sep- 
aration can be limited: pairs with separation > d cr it at two 
consecutive snapshots, but with separation < d cr it between 
these two snapshots, are automatically missed by the se- 
lection process. However, these instances are not expected 
to happen too frequently, and such pairs are not likely to 
penetrate close enough towards each other and have enough 
time to produce significant tidal damage. In any case, im- 
provements like these can be incorporated in future, more 
detailed analyses. 

It is interesting to concentrate also on the fraction of 
super-critical fate pairs relative to the sample of never- 
merging pairs (as opposed to the general population). The 
fractions are ~40%, 80%, 35%, 78%. This finding, along with 
the behaviour of the total-mass and mass-ratio distributions 
at the fate stage, support the idea that reducing d cr it aug- 
ments the influence of massive neighbours on pairs that 
never merge. The importance of super-critical pairs with 
small d cr it is more considerable at the fate stage than at 
the final stage. This is because the latter cases are a special 
subset of the former: they are the ones that are split apart 
more powerfully. 

This analysis suggest that single massive haloes (per- 
haps the main branch, but other massive secondary branches 
as well) enhance separations at the final stage and ultimately 
help split close pairs. The choice of adopted d cr it is very im- 
portant in this regard: small values of critical distance give 
more weight to the influence of these third-party massive ob- 
jects. At this point I speculate that further reducing d cr it will 
certainly mitigate the fraction of never-merging pairs, but 
those that never merge will be separated at even more vio- 
lent levels. On the other hand, while it is true that increasing 
Af cr it enhances mergers, the fraction of super-critical fate 
pairs is left largely untouched by this modification because 
there are not large numbers of supermassive neighbouring 
haloes to rip them apart at these impressive levels. 

4.2.4 The Lifetimes 

The lifetime of a halo pair goes from the initial to tf a te- ln 
general, the total duration of the pairs is very long, with 
over ~97% of the pairs lasting more than 1 Gyr. Notice 
that increasing M cr it increases the fraction of pairs that live 
more that 5 Gyr (from ~86% -> ~91% for the Close Sets, 
and from ~78% — > ~89% for the Very Close Sets - left versus 
right-hand panels in Figure [7|. Reducing d C rit, on the other 
hand, makes the total lifetimes shorter (upper versus lower 
panels in the Figure). It is possible that long total times are 
mitigated because adopting smaller critical distances tends 
to pre-select pairs that have a higher probability of merging. 

The approach period, which goes from imitiai to t cntTy , 
takes up a dominant fraction of the total time. In many 
ways, fapproach exhibits similar trends to itotai- Over ~95% 
of pairs have approach times greater than 1 Gyr. Increas- 
ing M cr it increases the fraction of pairs with approach time 
greater than 5 Gyr (from ~77% — ► ~84% for the Close Sets, 
and from ~66% to ~80% for the Very Close Sets). As is 
the case for the total time, reducing d C rit also shrinks the 
approach time. However, mergers do not happen at the end 
of the approach period. Nevertheless, adopting a small crit- 



ical distance selects pairs of haloes that are more prone to 
merge eventually, even if this happens long after the end of 
the approaching phase. 

As mentioned before, the total and approach times de- 
pend largely on the capabilities of the simulation. A simu- 
lation with either a larger box or higher resolution would in 
principle be capable of having even longer times (limited, of 
course, by the age of the Universe). A larger box allows more 
massive haloes with massive branches that are born earlier. 
Despite these limitations, it is important to report the ca- 
pabilities of the Millennium Simulation in regards to the 
evolution of close halo pairs. Higher resolution also allows 
all haloes to be born earlier, which can in turn be tracked 
down for longer periods, and more accurately. Even at late 
stages, small dissolving haloes in dense environments could 
be followed down for longer times (most likely making the in- 
teraction times l onger) if a simulation with h igher resolution 
were used (e.g.. iBovlan-Kolchin et all 12009). Alternatively, 
one can al so apply dynamical-time recipes to track these 
structures (jBinnev fc Tremaindll987h . 

Simulations also limit the short-time regime. In par- 
ticular, given that itotai is defined as the difference of 
two epochs associated with two simulation snapshots (i.e. 
itotai = tfate — ^initial) results in an artificial lower bound 
governed by the largest time snapshot of the simulation (for 
the Millennium Run, this corresponds to At = 0.0122 Gyr). 

The interaction time goes from the entry-time up until 
the last time the pair is within d C rit- For merging and sub- 
critical fate pairs, this ends at tfate ; otherwise it ends be- 
fore. Unlike the total and approach times, the distribution 
of interaction times has a noticeable short-time tail. This 
tail is more evident for the smaller choice of d cr it, where 
the prevalence of mergers suppresses the duration of these 
phases. For instance, the fraction of pairs with interaction 
times shorter than 1 Gyr increases from ~40% — > ~60% 
for the non-Massive Sets, and from ~30% — > ~50% for the 
Massive Sets. 

Using a large M cr it inhibits the fraction of pairs that 
interact for long periods because these pairs are more likely 
to spend most of their lives in the approach phase (see Sec- 
tion • In ^ cases, the contribution of pairs with dura- 
tions longer than 8 Gyr is negligible. 

Strictly speaking, the time spent within d cr it may be 
shorter than the interaction reported in this paper because, 
in many cases, the secondary halo may oscillate out and 
back inside the threshold on multiple occasions. Such in- 
stances are probably more important for subhaloe s orbit- 
ing a primary massive (perhaps central) halo (e.g., iBenson 
2005 1: iKhochfar fc BurkertllioOrj : lAngulo et al.ll2009l : IWetzei 
201 If ) . This level of detail is currently beyond the scope of 
this work, and can always be incorporated in more refined 
versions of the catalogue. 

4.3 Relation to Other Theoretical Work 

The nature of the evolution of pairs of haloes (or galaxies) 
from a cosmologically point of view remains relatively unex- 
plored. This section is devoted to discussing those few works 
that have attempted to rise beyond pre-prepared numerical 
experiments of a single nearly-isolated pair of galaxies. Here 
I simply discuss how their methods and assumptions com- 
pare to mine. However, a direct comparison of results is not 
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possible at this point because of the different set of goals of 
tho se papers and t he pr esent work. 

iTissera et all (|2002h are among the first research groups 
to attempt to simulate interactions from a cosmological per- 
spective. Smooth particle hydrodynamical (SPH) simula- 
tions are employed to study starbursts in merging galaxies. 
Their results suggest a double-triggering scenario where the 
first starburst is provoked early on in the interaction. Unfor- 
tunately, to attain such detail, they use a very small volume, 
Lbox = 5 Mpc h^ 1 (with h = 0.5) on the side, a nd only 64 3 
particles. Follow-up work by iPerez et all (|2006l ) uses a vol- 
ume with Lbox = 10 Mpc ft -1 and ~ 80 3 particles. The 
scope of this approach can be thought to be 'in between' 
the single-pair simulati ons and the present work. 

The approach of ISinha fc Hollev-Bockelmannl (1201 ll ) 
is very similar to mine, in the sense that they also fo- 
cus on dark-matter on ly simulations (unlike ITissera et al.l 
2002?; IPerez et al.l l2006. who include gas dynamics). These 
authors propose a notion called the 'halo interaction net- 
work', and suggest that fiy-bys at high redshift (z > 10) 
are more frequent than mergers. However, their simulation 
has Lbox = 50 Mpc on the side, only one thousandth of 
the volume in the Millennium Simulation. Their box size is 
comparable to the initial separations of pairs presented here 
(Figure |3}. It is not clear how representative their particular 
patch (or the merger histories within) actually is; or if they 
can produce enough massive background haloes where the 
m ajority of the hal o pairs happen to exist. See Appendix A 
of IWetzel fc Whitel (|20ld ) for a discussion on the impact of 
using a small simulation box on the clustering of subhaloes. 

iKitzbichler fc Whitel (120081 ) also employ the Millen- 
nium Simulati on, along with the se mi-a nalytic galaxy forma- 
tion m odel of ICroton et al.l (|2006l ) and iDe Lucia fc Blaizotl 
(2007). Their mock catalogue is extracted from a narrow 
light-cone with a field of view of 10 x 1.4 deg 2 . Their 
goal is very specific: to calibrate the time-scales connect- 
ing close galaxy pairs and their subsequent mergers. Given 
the nature of their objective, they are required to refine the 
stages immediately prior to the mergers. To do this, they in- 
clude additional dynamical friction times in their estim ates 
l|Binnev fc Tremaindll987l ; iBovlan-Kolchin et al.ll2008l ). Im- 
provements like this can certainly be implemented in the 
estimates shown in Fi gure [7) 

On a similar vein, iBerrier et al l (2006) use a simulation 
with Lbox = 120 Mpc h~ x on the side. Following halo merg- 
ers, these authors em ploy the analytic substructure model of 
IZentner et al. I (|2005h to track down subhaloes as they sink 
in their host groups. Their goal is also very concrete: to 
evaluate the number of close companions per galaxy. They 
argue that instead of using close halo pairs to predict the 
halo merger rate, the number of companions is a more suit- 
able quantity to test various galaxy formation scenarios. 
This research group then takes advantage of t hat approach 
to se arch simulated galaxy-pairs in isolation (|Barton et al.l 
l2007h and investigate star formation enhancements provoked 
by interactions. The types of selection criteria applied in 
these papers will be very valuable in future studies of halo 
pairs in different environments, and in more realistic appli- 
catio ns of my cat a logue . 

I Wetzel et al.l (|2008l ) also investigate how reliable close 
halo pairs are as proxies of mergers. However, they only 
concentrate on mergers of cluster-sized haloes (with masses 



> 5 x 10 13 M Q /i _1 ) at low redshifts (z < 1). To do this, they 
track pairs of 'halo objects' (i.e., the equivalent of an isolated 
halo in this paper), at the expense of ignoring subhaloes alto- 
gether. In a way, that work is complementary to mine: they 
manage to relate pairs that end up in different final groups 
but dismiss pairs of subhaloes within a given group. Pairing 
up all possible haloes and subhaloes at a given distance in a 
simulation, and tracking their evolution to test whether they 
merge or not, is a very computationally cumbersome prob- 
lem. A possible way to extend my approach in that direction 
is to compare branches that end up within a small bundle of 
groups in the local universe. However, at this point it is not 
clear how large the volume containing these groups should 
be. These difficulties illustrate how challenging the general 
study of close pairs as proxies of mergers can ac t ually be. 

On a slightly different subject. IWetzel et al.l (2008) also 
find that while rich environments have more halo pairs, 
mergers of such pairs are more likely to take place in less- 
dense environment. With the analyses presented in this work 
I argue that mergers are inhibited when the pairs involved 
happen to live in crowded environments. While it is true 
that this is reassuring, a direct comparison between the two 
works is impossible at the moment because their environ- 
mental scale is ~ 17 Mpc h~ l , while mine is restricted to 
regions that end up as single host haloes. 



5 SUMMARY 

This paper describes a statistical analysis of a very large 
catalogue of co-evolving halo pairs I constructed from the 
merger history trees in the publicly avai lable Millennium 
Simulation Database (|Springel et al.ll2005l ). 

The catalogue itself is comprised of over half a million 
pairs, selected by whether they come within a given crit- 
ical separation d cr it (Figure [T]). Only haloes that attain a 
mass larger than M cr i t at some point in their history are 
chosen (Figure [2]). Sub-catalogues are constructed to study 
the effects of modifying these criteria. The adopted values 
are d C rit = 1 and 0.2 Mpc h~ x ; and M crit = 8.6 x 10 10 and 
lO n M0/i -1 (corresponding to 100 and 1,000 simulation par- 
ticles respectively). See Table [1] for a summary of four sets 
arising from these combinations. 

In order to shed light on the evolution of these pairs, 
I establish a simple picture consisting of five basic stages: 
initial — > entry — » closest — > final — > fate (see Section [2.31) . 
Also, the following relevant time-scales are considered: the 
total, approach and interaction times (see Section [33} . 

Below I highlight some of the most interesting points 
we learn from this picture: 

• The majority of pairs begin at z > 4, with separations 
~ 10 Mpc h~ . The initial distance distribution is largely 
insensitive to d cr it, except for the subcritical pairs (with 
^initial < dcrit, only ~5% of them). 

• The distribution of the sum of the initial masses peaks 
at a few xlO lo M0/i _1 , regardless of ALrit- Initially, the two 
haloes in a given pair tend to be similar, with mass ratios 
typically greater than 1/10. 

• The small portion of sub-critical initial pairs is likely to 
correspond to the surviving sub-structures in cluster-sized 
host haloes. This is supported by the end-of-the-tail effects 
at low redshifts, high total-masses and low mass-ratios. 
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• As the evolution advances (entry — > closest — > final), to- 
tal masses increase steadily if they are > 1O 13 M0/i _1 , other- 
wise, they tend to decrease. This particular mass scale marks 
the range where tidal stripping becomes relevant. Mass ra- 
tios tend to decrease, making the discrepancy between the 
two haloes in the pair more evident. 

• Adopting a larger M cr i t shifts the mass distribution 
in the interacting stages to higher values: from 10 11 — > 
W 12 MqIi -1 . This also suppresses the high-redshift regime 
(with z > 4) significantly (for the stages after the entry 
time). 

• The final distance distribution peaks at d CI i t . Around 
40 — 50% of the pairs end up with super-critical separations 
(with dfinai > dcrit). Some of these pairs may be on their 
way to becoming permanently separated, while others may 
be temporarily on elongated orbits, bound to return together 
and merge eventually. 

• ~53% of the pairs in the primary catalogue never merge 
(these include pairs in which one of the haloes is absorbed 
by a third halo or pairs that do not merge by z = 0, but 
may merge in the future) . Reducing d cr it turns this number 
into ~42%; while enhancing Merit turns it into ~39%. With 
both effects, one gets ~27%. Over half of these non-merging 
pairs take place at z < 0.5. 

In general, merging pairs are more likely to be selected 
when the critical distance threshold is reduced. However, the 
relative fraction of never-merging pairs with super-critical 
fate separations is enhanced with this reduction. The same 
happens with the super-critical final separations. Squeezing 
the proximity criterion facilitates the selection of mergers, 
but also selects never-merging pairs that get ripped apart 
more violently. This hypothesis is supported by the fate 
total-masses and mass-ratios, which suggest a picture in 
which massive third haloes in the vicinity are responsible 
for preventing mergers. Indeed, it would be interesting to 
see if this effect remains with a simulation capable of follow- 
ing subhaloes with much higher mass resolution. Overall, a 
higher Merit also induces mergers. This is because super- 
massive third-party haloes in the vicinity capable of split- 
ting the pairs are less common. However, the violent super- 
critical effects are less important with this choice. 

The second aim of this paper is to provide rough esti- 
mates on the typical lifetimes of the pairs, and the duration 
of their interactions. Some of the main findings are: 

• On average, halo pairs last a few Gyr, spending most 
of this time in their approach phase. Interactions (from the 
first to the last time the haloes are within d C rit) last from a 
few hundred Myr to a few Gyr. Interactions are never longer 
than ~8 Gyr. 

• Adopting a smaller d cr i t promotes the fraction of merg- 
ers, shortening the typical interaction time as a consequence. 
Approach and total times are also reduced, but this effect is 
less noticeable. 

• Adopting a larger M cr i t produces longer lifetimes. This 
is because massive haloes are usually born earlier. However, 
this increase in duration happens mostly in the approaching 
stages (between the initial and entry times). 

In the discussion I mention that the nature of the initial 
stage depends strongly on the resolution of the simulation 
(which in turn affects the statistics of the approach and to- 



tal times). Nevertheless, it is useful to have an idea of the 
capabilities of the Millennium Simulation in regards to close 
pairs at the highest-possible redshifts. The mass-ratio distri- 
bution is also susceptible to the low-mass cut in the simula- 
tion. However, the types of applications the catalogue in this 
paper is ultimately intended for will also include necessary 
cuts in the types of haloes used. Lastly, resolution also affects 
subhaloes in the late stages of interaction. Higher resolution 
or the adoption of dynamical-friction recipes should, in prin- 
ciple, provide better estimates in this regime. In particular, 
interaction times are expected to be longer than those re- 
ported here. 

This paper marks the foundation of a series of papers 
dedicated to the analysis of this huge halo-pair catalogue. 
For the sake of space, several studies must be postponed for 
forthcoming papers. For instance, the basic quantities dis- 
cussed here (separation, redshift, total mass and mass ratio) 
constitute a four-dimensional parameter space, whose explo- 
ration is prohibitively unmanageable to be addressed in this 
exploratory paper. The same is true for the environment or 
kinematic quantities like relative velocities and spins. These 
studies are reserved for forthcoming papers. 

Ultimately, this catalogue will be the basis of future 
investigations of galactic interactions and double (and mul- 
tiple) quasar systems from a cosmological perspective. Such 
studies will require more detailed analyses, including the in- 
corporation of semi-analytic recipes and the handling of se- 
lection effects. For this reason, preliminary studies like the 
present work are a vital step. The results in this introductory 
paper will be essential pieces to guide this effort. 
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